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Abstract 

Transcription is regulated through interplay between 
transcription factors, an RNA polymerase (RNAP), 
and a promoter. Even for a simple repressive tran- 
scription factor that disturbs promoter activity at the 
initial binding of RNAP, its repression level is not de- 
termined solely by the dissociation constant of tran- 
scription factor but is sensitive to the time scales of 
processes in RNAP. We first analyse the promoter ac- 
tivity under strong repression by a slow binding re- 
pressor, in which case transcriptions occur in a burst, 
followed by a long quiescent period while a repressor 
binds to the operator; the number of transcriptions, the 
bursting and the quiescent times are estimated by re- 
action rates. We then examine interference effect from 
an opposing promoter, using the correlation function 
of transcription initiations for a single promoter. The 
interference is shown to de-repress the promoter be- 
cause RNAP's from the opposing promoter most likely 
encounter the repressor and remove it in case of strong 
repression. This dc-rcpression mechanism should be es- 
pecially prominent for the promoters that facilitate fast 
formation of open complex with the repressor whose 
binding rate is slower than 1/sec. Finally, we dis- 
cuss possibility of this mechanism for high activity of 
promoter PR in the hyp-mutant of lambda phage. 

Key words: transcription regulation; transcription fac- 
tor; transcription burst; transcription interference; math- 
ematical modeling 

Introduction 

The regulation of the activity of a particular gene in- 
volves a complex interplay between a promoter, an 
RNA polymerase (RNAP), and one or several tran- 
scription factors (TF) P, Ignoring the internal 
dynamics associated with transcription initiation, the 
probability for obtaining a successful RNAP elonga- 
tion initiation can be estimated from an equilibrium 



unbinding ratio of Tf[3, i]. When internal steps in 
transcription initiations becomes sizeable we need to 
consider the race between these steps and the kinetics 
of TF binding. 

The binding/unbinding rates of TF to bind to an 
operator is critically influenced by competitive non- 
specific bindings 0, 0] . Recent measurements of in vivo 
dynamics in an E. coli cell finds that a single lac re- 
pressor needs between 60 and 360 seconds to locate its 
operator Q. For TF whose copy number is of the or- 
der of 10 to 100 per cell, a cleared operator can remain 
free for up to about 30 seconds. In comparison, RNAP 
transcription initiation rates varies considerably, and 
can be as fast as 1.8 transcription initiation per second 
for a certain ribosomal promoter Q. Therefore, there 
is "room" for effects associated to the race between 
first bindings of a TF or an RNAP once the promoter 
is cleared. 

In a number of both procaryotic and eucaryotic sys- 
tems, the promoter activity are not only influenced by 
TF, but are also modulated by interfering promoters 
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[12, lia, llj, |15|. For example, the reg- 
ulation between lytic and lysogenic maintenance pro- 
moters in the P2 class of bacteriophages involves tran- 
scription interferences (TI) as well as TF's that repress 
the promoter activities [ll|. And in lambdoid phages 
the initial lysis-lysogeny decision is modulated by TI 
between the promoter PRE activated by CII and the 
promoter PR repressed by CI. 

Dodd et al. [ISj] presented a framework to deal with 
TI and multiple TF's, using an assumption about fast 
equilibrium reactions of TF-binding and closed com- 
plex formation. In the present paper, we develop a for- 
malism that deals with the competition between time 
scales of TF binding/unbinding and transcription ini- 
tiation process, and examine the effect of interference. 

FiglD shows a single promoter pS with an oper- 
ator site for a repressive TF (left panel), and with 
a convergent promoter pA (right panel). For both 
cases, we illustrate the three basic steps of transcrip- 
tion initiation: (i) RNAP reversible binding to form 
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Single Step Model 



Figure 1: (a) Model of promoter pS with a single TF that 
represses the promoter by competitive binding to an opera- 
tor that overlaps with the promoter. The promoter activity 
is given in terms of the three-step Hawley-McClure model 
for transcription initiation, with indicated transition rates 
for formation of closed complex, that of open complex, and 
elongation, (b) Same as in (a), but with addition of a con- 
vergent promoter that interfere with both RNAP binding 
to pS and with binding of TF. 



a closed complex, (ii) irreversible transition to open 
complex, and (iii) initiation of transcription elongation. 
The rates for these three steps are promoter depen- 

3, 



dent 
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As for the initial binding, given the 
fact that the maximum activity for ribosomal promot- 
ers reaches 1.8 transcriptions per second Q, tire time 
needed for an RNAP to diffuse to a promoter cannot 
be longer than ~ 0.5 sec. Regarding the later steps 
where RNAP forms open complex and subsequently 
initiates transcription to leave the promoter, their time 
scales may vary a great deal from one promoter to 
another [li, ^ ^M,^- 

In the following, we will investigate in detail how 
these time scales play together to determine the extent 
to which a promoter is sensitive to repressors and to 
clearance due to the interference by elongating RNAP's 



from other promoters 24 



Models 

We study the promoter activity under influence of tran- 
scription factor(TF) and transcription intcrfcrence(TI) 
based on mathematical analysis on simple models of 
promoter in tire following three levels. Our goal is 
to understand regulation of the three step model for 
transcription initiation originally proposed by Hawley- 
McClure[ll,[l3, but we also analyze its simplified ver- 
sions, i.e. the single step model and the two step model. 
The comparison of these three levels of models gives us 
intuitive understanding of the promoter behavior. 



Two Step Model 



Three Step Model 



Figure 2: Schematic illustrations for the single step model, 
the two step model, and the three step model of the elon- 
gation initiation. 

Three Models for Elongation Initiation 

Let us start by describing the bare models with neither 
TF nor TI (FigO. 

i) The single step model of transcription initiation 
is the model where the whole process is dominated by 
a slowest step, thus its elongation initiation is repre- 
sented by a simple Poissonian process with the rate 

no. 

ii) In the two step model, the transcription initia- 
tion consists of two steps: first, RNAP binds to the 
promoter site with the on-rate kon, and then initiates 
elongation with the rate ke- The transcription initia- 
tion rate for the overall process flo is given bv[l3| 



no - 



h h 



1 



(1) 



with 



(2) 
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'''on = T ; = T~ • 

/Con "-e 

The last expression of ([T|) simply shows that the av- 
erage interval of elongation initiation l/^o is the sum 
of the two times: Ton, the time for RNAP to form the 
on-state, and r^, the time to start elongation in the 
on-state. 

iii) In the three step model, two states within tire 
RNAP binding state arc differentiated: the one with 
closed DNA complex and the other with open DNA 
complex. The transition between the RNAP unbind- 
ing state (off-state) and the RNAP binding state with 
closed DNA is reversible, and characterized by the bind- 
ing rate fcf, and the unbinding rate When RNAP is 
in the closed complex state, the transition to the open 
state is irreversible with the rate ko. Finally, the open 
complex is followed by elongation initiation with the 
rate fcg. This three step model of transcription initi- 
ation was originally proposed by Hawlev-McClureflil 

The three step model reduces to the two step model 



with the effective on-ratc k*^^ given by 

k* EE ^° 

l + ku/kb 



(3) 
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in the case where the off-state and the closed DNA 
binding state are in equiUbrium. This is fulfilled when 
the initial reversible process of RNAP binding/unbinding 
is faster than the other processes: fcf,, fc„ 3> fco, [l^ . 
The effective on-rate in eq.Q can be understood 
as the open rate ko reduced by the equilibrium expec- 
tation of being unbound. 

The overall elongation rate Hq for the three step 
model has been shown [l4| to be 



(a) Bare promoter activity 



1 



1 



1/kb + 1//C*„ + l/ke 



n 



with 



n 



1 

kb ' 



ko 



1 

ho 



(4) 



(5) 



The time r* is the time for the system to form an open 
complex after an RNAP binds to form a closed state for 
the first time. It is the sum of the two times: (i) Tq, the 
time to form an open complex without unbinding, and 
(ii) the binding time rt multiplied by the average num- 
ber of times of RNAP unbindings before forming an 
open complex, ku/ko (a detailed explanation of math- 
ematical interpretation is given in the appendix of the 
supplement). Note that this expression holds for a gen- 
eral case, not limited to the case where the two step 
approximation is valid. 

In the above discussion, we have ignored the self- 
occlusion effect, where the next RNAP cannot bind to 
the operator site until the previous RNAP goes away 
from it. If we include this self-occlusion effect, the bare 
activity flso should be 



(6) 



with Tso being the time that RNAP needs to clear the 
promoter. 

Transcription Factor 

For each of these models, we consider the effect of a 
repressive transcription factor (TF), which we assume 
completely prevents RNAP from binding while it binds 
to the operator site. It is also assumed that RNAP 
binding to the promoter site prevents TF from binding 
to the operator site. The binding and unbinding rates 
of TF are denoted by k"^^ and fc^^ ' respectively. 

We will study, in particular, the strong repression 
regime, i.e. the dissociation ratio k'^^ /kj^ is small. In 
such a case, TF binds for most of the time, prevent- 
ing transcription initiation, but once a TF falls off, 
the promoter is free to initiate a burst of transcription 
elongations until another TF binds to the operator site 
(FigHb)). 



(b) Promoter activity with TF 



(C) Promoter activity with TF & TI 
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Figure 3: Schematic diagrams for the time sequence of 
a promoter activity for a bare promoter (a) , a promoter 
with TF regulation (b), and a promoter with TF un- 
der TI (c) . The vertical lines represent the times when 
transcriptions are initiated. The shaded(cyan) inter- 
vals labeled as TF represent the time intervals when 
a TF bounds to the operator site, thus the promoter 
cannot initiate transcription. Under the TF regulation, 
the transcription bursts take place while a TF does not 
bind. The arrows indicate the times when interfering 
RNAP's from pA arrive at pS and remove both TF and 
RNAP at pS; TI triggers transcription bursts. 

Transcription Interference 

The effect of transcription interference(TI) on the pro- 
moter pS is examined by exposing it to transcribing 
RNAP's from another promoter pA in parallel Q or 
in convergent 0, |Tl| configuration (the latter case is 
illustrated in FiglUb)). The interfering promoter pA 
is characterized by the transcription initiation rate 
and the initiation interval distribution pa{t). The RNAP 's 
from pA are assumed to clear both the promoter and 
the operator sites of pS (sitting duck interference) and 
to occlude them while passing. This causes bursts of 
transcriptions after the interference until another TF 
binds(Figl3i;c)). 

There are several additional complications related 
to TI. (i) The RNAP sitting at the operator and the 
TF at the promoter of pS may not simply fall off by 
the interfering RNAP from pA, but may block it (road- 
block effect), (ii) Between the promoter and the opera- 
tor, there should be time difference for the sitting duck 
interference and the occlusion to take place because 
they extend over a certain finite size and are located 
at difference places along DNA. (iii) The interference 
may also take place through collision with an RNAP 
from pA after an RNAP from pS starts elongation, (iv) 
The interference between pS and pA should be mutual, 
namely, pS can also interfere in the pA activity while 
pA interferes with pS. 

In the case where pS and pA arc in a parallel config- 
uration, the collision effect (iii) and mutual interference 
(iv) do not exist. Even in a converging configuration, 
the collision effect is not significant when the distance 
between pS and pA is short, i.e. the traveling time 
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between the two promoters is much shorter than the 
activity interval of the promoters. As for the mutual 
interference, the effect of pS on pA is negligible when 
the activity ^Ia of pA is much larger than the activity 

n of pS. 

These effects (i)~(iv) introduce further complica- 
tions in the problem, but we are going to ignore all of 
them in the following. 

Outline of Theory 

The quantity wc arc going to examine is the averaged 
elongation initiation rate, or activity of pS, under the 
influence of TF and TI. Under the repression by TF, 
a promoter initiates transcriptions in bursts and we 
will see how TI can activate the promoter. This effect 
can be prominent especially when the TF repression 
is strong and the time scale for TF is slow. In this 
section, we outline the theory. Detailed derivations of 
formulas are given in the supplement. 

Single Promoter Property 

As tools for the analysis, we use the following two func- 
tions: (i) p{t) , the probability distribution for time in- 
tervals between subsequent elongation initiation events, 
and (ii) C(i), the averaged time-dependent rate of elon- 
gation initiation after both the promoter and the oper- 
ator sites are cleared. We first examine p(r) and C{t) 
for pS without TI, but under the effect of TF. 

The average elongation rate fl without TI is the 
inverse of the average elongation interval, thus it is 
related with p(t) as 



n = 



p{t)t cLt 



(7) 



The time dependent elongation rate C(t) is actually 
a correlation function of elongation initiations without 
TI because it can be regarded as a probability density 
of initiation at the time t provided that there was an 
initiation at t = 0. This can be directly calculated 
from p{t). For large i, C{t) approaches the promoter 
strength fi, 

n = lim C{t) (8) 

t — >oo 

because the effect of the initiation at < = lasts only a 
finite time. 

Transcription Interference 

Now, we consider TI. Under the influence of interfering 
promoter pA, the promoter pS and its operator site 
are assumed to be cleared every time an RNAP from 
pA passes, and the activity of pS will change as C{t) 



after that. Thus the time averaged activity during the 
interval of length r is given by 



If 

T- Jo 



C{t) dt, 



(9) 



where we have included the occlusion time Tqcc- The 
occlusion time Tqcc (= 1 2 sec) [25| is the time where 
the pS promoter cannot bind a new RNAP due to 
a transcribing RNAP from pA. This effect is not in- 
cluded in the correlation function C(t), because the 
correlation function defined here is a single promoter 
property. 

The overall average activity of pS is the average 
of eq.® over the interval distribution of pA, namely, 
Pa{t). It is important to notice, however, that this 
average is not with the weight pAi^) itself but with the 
weight proportional to pa{t)t because the probability 
that a given time falls in the interval of length t is 
proportional to pa{t)t, not pa{t). Therefore, the final 
expression for the elongation rate under TI is 



/ Pa{t)t 




— Tocc 

C{t) dt 


dr 


J Tocc 


T Jo 




r 

Jo 


Pa{t) 


T dr 





(10) 



The occlusion effect by RNAP from pA is explicitly 
included as a finite Tqcc, but the self-occlusion effect, 
that the RNAP from pS blocks its own promoter site 
pS, should be included in the correlation function C{t) 
if it is considered. 

In addition to ignoring (i) roadblock effect, (n) time 
difference between the promoter and the operator, (iii) 
RNAP collision, and (iv) mutual interference, we will 
further approximate pA as Poissonian, namely. 



PAir) = r^Ae 



(11) 



and also ignore the occlusion time by putting Tqcc = 0, 
and the self-occlusion effects. 

Evaluation of p(r) and C{t) 

By assuming each elementary process, such as binding, 
unbinding, elongation, etc., to be a Poissonian process 
with a given rate, we can obtain analytic expressions 
for p(t) and C (t) , from which we can calculate the over- 
all elongation rate f2 for pS for various situation with- 
out TI. Using these functions, the elongation activity 
under the influence of TI is estimated from eq. pUj) with 

'^occ 0. 

Detailed derivation of mathematical formulas is given 
in the supplement. In the following, we will describe 
results obtained from those analytic expressions. 
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(a) Single Step 




(c) Three Step 
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Figure 4: The interval distribution p(r) (dashed green 
hnes) and the time dependent elongation rate C(t) (solid 
red lines) for the single step (a), the two step (b), and 
the three step model (c). The parameters for the three 
step model are ki, = ls~^, ku = ls~^, ko — 0.2s^^, 
and fee = 0.1 s~^, which gives the average elongation rate 
fio = (l/feb + 1/feon + 1/fee)"^ = 1/21 s-1 and the effec- 
tive on-rate feon = feo/(l + feu/fefc) = 0.1 s~^. The param- 
eters for the two step model are determined so that they 
behave similarly, i.e., the on-rate feon = feon a-nd the elonga- 

(2) 

tion rate in the on-state fee for the two state model given 
by ki^^ — [l/flo — 1/feon) The insets show the behav- 
iors around t ~ with the asymptotic curves (dashed blue 
lines) . 



Results 

We present the numerical evaluations of our expres- 
sions for various situations to clarify dynamical effects 
of TF and TI on the promoter activity. 

Activity of a Bare Promoter 

Let us start by comparing the three models in a bare 
form, i.e. without TF and Tl. 

FigH] shows the elongation initiation interval dis- 
tribution p{t) (dashed green lines) and the time depen- 
dent activity C{t) after the promoter site have been 
cleared by the competing activities (solid red lines). The 
parameters are chosen for the three step model, and 
those for the two step and the single step models are 
determined to match them with the three step model 
using eqs.(l3|) and ([4]), namely, fcon = fcon ^^'^ to give 
the same overall activity fio for all the cases. 

In the single step model, the elongation initiation is 
Poissonian, and the interval distribution p(r) is a sim- 
ple exponential with the elongation rate fio- As there 
will be no correlations between subsequent initiations, 
the activity C{t) is given by the constant fio- 

In the two step model, p{t) and C{t) rise linearly 
from zero as konket (the inset in FigHJ^b)). This is 
because the RNAP has to bind to the promoter site 
with the rate fcon before it initiates elongation with the 
rate fcg- The difference from the single step model is 
seen in the time scale t < min(l/fcon^ 1/fce)- The two 
step model reduces to the single step model in the case 
either fcg 4C fcon or ke 3> fcon, but these two cases show 



quite different behaviors in reaction to TF, as we can 
see in the following subsections. 

In the three step model, the promoter goes through 
two states after RNAP binding. Therefore p{t) and 
C{t) increases initially as {\/2)ki,kokJ? around t = 
(the inset in FigHI^c)). In the case of fast equilibration 
in the initial transition(fcti, fc^ ^ fco), the three step 
model reduces to the two step model with an effective 
on-rate k*^ given by eq. ([3]) . 

In general, the main feature of an increased num- 
ber of intermediate RNAP-promoter states causes an 
initial rise of p{t) and consequently C{t) to be of in- 
creasing order in t or t. Also the peak in p{t) becomes 
sharper, which in principle could give a non-monotonic 
behavior of C{t). For any realistic parameters, how- 
ever, we find monotonic C{t) for the promoters without 
TF. 

Activity of Regulated Promoter by TF 

We now consider a promoter which is regulated by a TF 
that acts as repressor as illustrated in Fig[T]Ja). Under 
strong repression by a slow binding TF, transcriptions 
occur in bursts with quiescent periods of the length 



ttf 



1 

Ttf ' 



(12) 



when a TF binds to the operator and suppresses the ac- 
tivity. We will see the general expressions of promoter 
activity IItf repressed by TF can be put in the form 
that allows direct interpretation in terms of transcrip- 
tion burst. We evaluate the time-dependent activity 
C{t) for various parameters under the influence of TF, 
whose binding and unbinding rates are fc^^ = Is""'^ 



with kl^ Ik'^^' = 0.1, 0.01, and 0.001. C(t)'s without 
TF and with TF which never unbinds, i.e. fc^^ = 0, 
are also plotted for comparison (dashed green lines). 

Single step model 

The TF effect on the single step model is rather straight- 
forward. The expression for C{t) is given by 



C{t) 



no 



, TF 



kr 



jLTF 



I TF - 

h e 



(13) 



which is plotted in FiglSja). Immediately after the 
promoter is cleared at < = 0, the promoter activity re- 
covers to the bare value f2o: but the initial high activity 
decreases as a TF binds around t ^ 1/kJ^. In the lat- 
ter stage, the transcription initiation is determined by 
the equilibrium probability of having a free promoter, 
/i^^I^ + ^u^)- Therefore, C{t) shows an exponen- 
tial decrease from the initial bare activity fip to the 
repressed level of averaged activity. 



k 



TF 



TF 



fcTF + fcTF 



(14) 
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- (a) Single Step 













10' 10- 10^ 
1 [sec] 




10^ 10- 
[sec] 



Figure 5: The time dependent activity profile C{t) with 
TF for the single step (a) and the two step models with 
kon = 0.1s~^(b) and ls~^ (c). The bare activity is one 
transcription initiation per 20 sec: fio — 1/20 s~^. For 
each case, we show five curves: the un-repressed case with- 
out TF (top dashed green lines) , the repressed cases by TF 
with the binding rate k^^ = ls~^ and the unbinding rate 
fcjf^ — 0.1 s~^ (top solid red lines), 0.01 s~^ (middle solid 
red lines), 0.001 s~^ (bottom solid red lines), and with TF 
that never unbinds (bottom dashed green lines). The ar- 
rows indicate Cmax given by eq. (|22p and tpi given by ea. (|23|) 
for the two step models. For the single step model(a), the 
RNAP activity is limited only by a binding event once ev- 
ery 20 sec. For the case (b) of the two step model, each 
step takes 10 sec, while in the case of (c) where the on-rate 
is fast, overall activity limited by an elongation initiation 
time of 19 sec. 

for t ^ Note that this simple "equilibrium re- 

pression formula" Q for transcription repression holds 
only for the single step model. More subtle competi- 
tion comes into the problem for the two and three step 
model, as we will see below. 

It is interesting to see that the equilibrium formula 
can be also put in the form 



r^TF — 



nhst 



with 



Tbst 



Tbst + TTF 



(15) 



(16) 



and ttf defined in (fT2|) . This allows direct interpreta- 
tion in terms of transcription burst; Tbst and ribst are 
the typical time scale and the number of transcriptions, 
respectively, of a single transcription burst, and ttf is 
the typical time scale of the quiescent period between 
the bursts with TF bound to the operator. The expres- 
sion (|15p represents that the average promoter activity 
Htf is given by the number of transcriptions in a burst 
nbst divided by the time interval between the consecu- 
tive bursts, Tbst + ■''TF- Note that the expression (fT5|) 
itself is valid in general case and not limited to the case 
where the transcriptions occurs in burst, namely, the 
promoter is strongly repressed by a slow binding TF. 

Two step model 

The situation is a little more complicated for the two 
step model. In FiglSI^b) and (c), two cases are shown: 



one with fcon = 0.1 s~^ and the other with fcon = 1 s~^; 
In the first case, the time scales of the two transitions, 
the on-rate and the elongation rate, are same, but in 
the second case, the on-rate is much faster than the 
elongation rate. The elongation rate fcg are chosen to 
give the same bare activity ilo for the two cases. 

The general behavior of C{t) is that (i) first it in- 
creases as konkgt until TF starts binding, (ii) then it 
reaches a plateau value, and (iii) finally it goes to the 
steady activity fixF averaged over long time. 

The time averaged activity with TF is given by 



UTF 



TF 



[fce/(fco 



(17) 



with r^o being the bare activity of the two step model 
([T]). Note that the repression factor by TF, i.e. VIty/^o, 
is given by the "equilibrium formula" (|14p only when 
fee 3> fcon- In the other limit, TF cannot repress the 
promoter as one might expect from the dissociation 
constant of TF, /k'^^ . 

This time averaged activity p7)) can be also ex- 
pressed in the same form with eq. ()15|l . 



"bst 



TF 



Tbst + Ttf 



but TT-bst and Tbst are given by 

k: 



n-bst 



^on 
lTF • 

^6 



Tbst 



jlTF 



n-hst - 



(18) 



(19) 



Here, ribst can be understood as the number of tran- 
scriptions in a burst before a TF binds to the operator 
because kon/kj^ is the "winning ratio" of RNAP to 
TF for binding. The bursting time Tbst is the sum 
of the binding time of TF, 1/kJ^, and the elongation 
time, 1/fce, multiplied by the number of transcriptions. 
Again, this expressions is valid in general case although 
it is interpreted best in the bursting situation. 

In the strong repression limit where the bursting 
time is negligible compared with the quiescent time, 
we have 



n 



TF 



^on 
TT 
b 



■k 



TF 



when Tbst < TTF. 



(20) 



Note that the time averaged activity in this limit does 
not depend on the elongation rate fcg in the on-state. 
This is because the time scale is set by the slowest 
rate k"^^ . The promoter produces a burst of ribst (= 
kon/kj^) transcriptions while a TF is not bound, but 
once a TF binds, it has to wait a time ~ ttf(= V^u^) 
for TF to unbind. 

In the case fc^^ ^ fce ^ fcw^i the plateau becomes 
a maximum; C{t) can be approximated as 



Cit) 



(21) 
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for i < l/Zce • ln(Ae/^«^) (see the supplement). From 
this expression, we can estimate the maximum value 

^on IT'hst 



Three Step Model 



Cn 



with the plateau time 
tp\ = 



nhst 



(22) 



(23) 




10" 10' 10- lo-" 
/ [sec] 



■(b) 








. / kt,=1 s''' 

/ kj,=10s''' 
kg=0.05s"'' 






10"' 10** lO' 


10 


' 10- 



t [sec] 




10 ' 10" 10' 10- lO-" 
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for the time region 



kr) 



<t <t 



pi- 



(24) 



The time dependent activity C(t) shows maximum 
after the promoter site is clarified. The maximum value 
(|22p can be understood as ke times branching proba- 
bility to the on-state fcon/(fc(f^ + fcon)- This and eq.([5D]) 
show that the promoter repression by TF is determined 
by the competition between TF and RNAP for bind- 
ing to DNA, namely, between the binding rate kj^ and 
kon- Therefore, even if the bare activity is the same, 
the repression by a TF can be quite different. This can 
be seen in FigO /con = 0.1 s~^ (b) and Is^^ (c) with 
the same fio = 0.05 s~^. The repression in FigHJc) is 
about 10 times weaker than that in (b), because fcon is 
10 times faster. 

After TF falls off from the operator site, the pro- 
moter produces a burst of nbst(= kon/kj^) transcrip- 
tions on average before another TF binds. Note that 
tpi « Tbst in the case fcon > fc^^, namely, Tibst > 1- 



Three step model 

In the full three step model, the RNAP have to pass 
through a closed DNA complex state first. The transi- 
tion between this closed complex state and the off-state 
is reversible, but its rates can be relatively fast com- 
pared with the transition rates of the following steps. 
The fast initial binding process tends to make TF re- 
pression less efficient. This has been verified by mea- 
surements on promoters with strong RNAP binding 
affinity 

The general expression for the time averaged activ- 
ity with TF is again given by 



r^TF — 



?lbst 



with 



"-bst 



Tbst 



Tbst + TTF 



ko 



kb 

fc() fco ~t- fc ji 



I,TF 



ribst I TT 



(25) 

(26) 
(27) 



The number of transcriptions in a burst ribst is now 
given by the winning ratio kb/kj^ of RNAP to TF 



Figure 6: The time dependent activity profile C{t) with 
TF for the three step models. For each case, we show 
five lines: the un-repressed case without TF (top dashed 
green lines), the repressed cases by TF with the binding 
rate fcjf = 1 s~^ and the unbinding rate fcjf^ = 0.1 s~^ (top 
solid red lines), 0.01 s^^ (middle solid red lines), 0.001 s^^ 
(bottom solid red lines), and with TF that never unbinds 
(bottom dashed green lines). The arrows indicate Cmax 
given by ea. (|34p and tpi given by ea. (|30p . For all cases, the 
bare activity is Qo ~ 0.05 s~^, but we focus on the promot- 
ers with fast open complex formation, namely, the larger 
effective on-rate k*^ = 10/11 s~^. The case (a) corresponds 
with a strong closed complex binding ku/kb 1, whereas 
(b) and (c) deals with a weakly binding RNAP. The differ- 
ence between (b) and (c) illustrates the effect of a 10 times 
faster RNAP binding rate to the promoter. 



multiplied by the branching ratio fco/(fco + ^u) in the 
closed state to the open state. The bursting time Tbst 
is the sum of the TF binding time 1 / kj^ and the time 
needed for elongation after RNAP binding to the pro- 
moter multiplied by the number of transcriptions. Note 
that the bare activity fJo in eq.(|4]) can be expressed as 



'tbst 
Tbst ' 



(28) 



which also holds for the other two models. 

It is easy to see from eqs.(j4|) and (|25l) that the re- 
pression factor r^TF/f^o is given by the "equilibrium 
formula" (fTil) only when k^ , fco ^ fcfc , kj^ , namely, the 
internal time scales are negligible. Note that the ex- 
pression (|25p can be put also in Michaelis-Menten form 
using (effective) dissociation constants (See Appendix). 

In the strong repression limit where the bursting 
time is negligible, it is easy to see that 



n 



TF 



TTF 



kn 



kb 



l,TF 

'^b 



■k 



TF 



ku ~t~ ko 
when Tbst < TTF 



(29) 



from eqs.(|25l) and ([26]). This reduces to eq.(|20l) with 
fcon replaced by k*^ of eq.([3]), in the case of a weakly 
bound closed complex (fc^ ^ kb, ko), because fcon ~ 
kokb/ku in this limit. 

FigjS] shows the time dependent activity profiles 
for three promoters whose bare activities are similar 
but with different closed complex formation transition 
rates. The first two cases, (a) and (b), are for the 
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same k*^ = 0.909 s~^, but for the last case (c) fc*,^ = 
0.0545 s~^. One see that the promoters respond dif- 
ferently to repression by a TF. The arrows show the 
maximum value Cmax of eq. (|22p with the plateau time 



and ribst for the three step model. 



(30) 



Schematic Description for Time-Dependent Ac- 
tivity 

With all these results, Figl7]summarizes the behavior of 
time dependent activity C{t) for the promoter with fast 
initial binding/unbinding under the strong but slow TF 
repression: 

(fc„, kh) > (fc^^, fco, ke) ^ ^u^i (31) 

where we have a typical bursting of transcriptions with 
n-hst > 1, ttf > V ~ Tbst- (32) 

After the clarification of the promoter and operator 
sites, the activity increases initially as 



C{t) w ^kbkoket'^, 



for t < -ip (33) 



until TF starts binding. 

Then, it reaches the (maximum) plateau value: 



nhst 



n-bst 



1 



pi 



"bst + 1 1/fco + l/fce 



for t < t 



pi- 



(34) 

Finally, C{t) diminishes down to the long time av- 
eraged steady value with TF, 



Tbst + ttf 



ttf 



(35) 



From eqs. ([M)) and (pS)) . the enhancement factor 
/cnh that the promoter can be activated after the clear- 
ance of the site is given by 



enh 



Cmax ^ T hst + ^TF ^ ^TF 



(36) 



This expression formalizes our original discussion that 
one obtain large relative peak activity when TF repres- 
sion is strong, (xbst, *pi) <C ttf, but slow k]^^ <^ kb, 
namely, n^st 3> 1. The promoters with shorter "inter- 
nal time" 1/ko + l/ke have larger relative peak activity, 
and therefore they will be more prone to de-repression 
by TI. 



C(t) in log scale 



Without TF 




1 n 



bst 



o u b 



Ipl 



t in log scale 



Figure 7: Schematic activity profile of a promoter, with 
and without a TF that acts as a repressor by occluding the 
promoter site. Without repressor the promoter activity is 
set by the time that the promoter takes to pass through 
the three steps to initiate elongation, whereas a repressor 
reduces the final promoter activity to the extent propor- 
tional to the dissociation ratio k^^ /k'[^ . Shortly after the 
promoter clearance, the activity recovers to reach the maoc- 
imum value ribst/ipi until t < ipi. This can be much higher 
than the steady activity nbst/TTP when tt\ <^ ttf, i.e. the 
promoter is strongly repressed (ttf ^ tpi,Tbst) by a slow 
binding TF{k'^^ < kb or ribst > 1). 

Interfering with Regulated Promoter Ac- 
tivity 

Wc now consider the interfering promoters pS and pA 
where pA is relatively strong in comparison with pS, 
and pS is strongly repressed by a slow TF. In this case, 
the average activity JIti is given by eg. pU]) . using the 
time-dependent activity C{t) without TI and the elon- 
gation interval distribution paIt) of pA. 

In the following, we ignore the occlusion time Tqcc 
by RNAP from pA; This should not be bad for the pro- 
moter whose activity is of order or less than 0.1 s^^, but 
may not be so good for a more active promoter. For 
Pa{t), we will use the exponential distribution pip , 
which corresponds to the single step Poissonian pro- 
moter pA. 

The expression for I^ti of pO|) basically gives the 
average of C(t) over the typical time scale of pA, which 
is fl^^. Therefore, if you look at fixF as a function of 
Qa, then Oxi would show a maximum around CI a 
l/^max in the case that C{t) has a maximum around 

t ~ tmax ■ 

Two step model: we can obtain the explicit expres- 
sion for Qti, which can be approximated as 



nhst 



TI 



l,TF 



tpi + n- 



I.TF 
H 



TIa 



(37) 



for Qa ^ in the regime kj^ 3> fee S> ^'u^■ Here, 
ribst is the number of transcriptions in a burst defined 
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in eq. p^ . and tp\ is the plateau time ([^5]) . This ex- 
presses the activity of the de-repressed two-step pro- 
moter in terms of the product of two factors; the first 
factor corresponds to the averaged activity for the burst- 
ing whose interval is given by tpi+Q~^^ . This factor rep- 
resents the de-repression by the interference through 
removal of TF by RNAP from pA before it dissociates 
by itself. The second factor represents the suppression 
by removing the open complex, i.e. sitting duck inter- 
ference. 

The expression p7p shows a maximum 



Three Step Model 

(a) k^^l.k^^lO. k^^l (b) k^^l, l^,^!, k^^O.2. k^^.l (c) k,,^!. k^^^lO, k^^]. k^.^.05 



rihst 



TI,max 



at flA 



Pl 



t 



pl 



(38) 



which corresponds to cq. (|22p . The promoter activity 
is actually increased by the transcription interference. 
The enhancement factor, or the ratio r^Ti.max/^^TF, is 
the same as in eg. ipS)) . 

Three step model: We cannot write down a com- 
pact expression, but FigIS] shows Q^i vs. ^Ia (lower 
panels) along with the corresponding C (t) (upper pan- 
els) for the three step model with different values of 
parameters. One can see the correspondence between 
the upper and lower panels: The activity as function 
of approximately resembles the activity profiles a 
time t ^ l/f^A plotted in corresponding upper panels. 
Notice also that the potential activation by a conver- 
gent promoter is largest for large and fcg, as ex- 
pected from ea. ([5S)) . Finally, the relative effect of de- 
repression can be very large, in the case of very slow 
dissociation rate for the transcription factor. 

Summary and Discussions 

We have presented a mathematical framework that ex- 
presses the dynamics of a promoter in the Hawley- 
McClure model. The formalism opens for discussion on 
the promoter activity with transcription factors (TF) 
and transcriptional interference (TI) by an interfer- 
ing promoter, and allows us to deal with the interplay 
among these elements. 

In our formalism, the activity of a single promoter 
is characterized by the correlation function C {t) , which 
represents the averaged time dependent activity after 
the transcription initiation at t ~ 0. Any modifica- 
tions to the activity of the single promoter, such as 
that by TF, are taken into account through the cor- 
relation function C{t). On the other hand, the ef- 
fects from an interfering promoter punctuates the pro- 
moter/operator activity with the time scale of tran- 
scription initiation from the interfering promoter. This 
is represented by the expression pO|) . 

We have studied the effects of a repressive TF on 
the promoter activity. The general expressions for the 




Figure 8: The time dependent activity profile C{t) with TF 
(upper plots) and the average activity Qti with TF and TI 
vs. Q,A (lower plots) for the three step models. For each 
case, we show five lines: the un-repressed case without TF 
(top dashed green lines), the repressed cases by TF with 
the binding rate kj^ = 1 s~^ and the unbinding rate kl^^ = 
0.1 s~^ (top solid red lines), 0.01 s~^ (middle solid red lines), 
0.001 s^^ (bottom solid red lines), and with TF that never 
unbinds (bottom dashed green lines). The arrows indicate 
Cmax given by ea. p4)l and tpi given by eg. pOt for the upper 
plots, or Slxi.max and Qa by ea. (|38fl with fcon replaced by 
fcon and ribst and tpi replaced by ones of the three step model 
for the lower plots. Notice that the activity at a given level 
of pA activity, reflects the average activity of the promoter 
up to a cut-ofi^ time of I/^Ia- 



promoter activity can be put in the form that is as- 
sociated with the transcription burst, namely, a burst 
of ribst transcriptions during the bursting time Tbst fol- 
lowed by a quiescent period of the length ttf- This 
is actually what happens in the case of a strongly re- 
pressed promoter by a slow binding TF. It should be 
noted that the "equilibrium formula" ([T^ for the pro- 
moter activity repressed by TF is not valid unless the 
time scales of internal processes are negligible com- 
pared with binding/unbinding times of RNAP, because 
a TF competes with RNAP for binding to DNA. 

Under the transcription interference (TI) consid- 
ered in the present work, an interfering RNAP simply 
clears both the promoter and operator sites. If the pro- 
moter is strongly repressed by a TF, such interference 
is most likely to relieve the promoter out of repression, 
and interrupts the quiescent period to shorten to 1/ Qa 
when ^Ia > k'^^ ■ 

Experimental Observations 

Let us discuss experimental relevance of our theoretical 
results for transcription burst and its modification by 
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transcription interference. 



Experimental Proposal 



Transcription Bursts: Experimentally, bunched pro- 
moter activities have been seen in several eucaryotic 
systems involving TF's (26l . 27 1, and they have been 



interpreted to occur in response to transformation in 
heterochromatin states, or as a result of a promoter 
approaching transcription factories [2^. For procary- 
otes, bunched activities have only been observed for the 
promoter -Piac/ara under the fully induced condition [2^, 
and has been interpreted without TF[3^. 

On the other hand, transcription bursts induced by 
activators have been examined in models and exper- 
iments on the yeast GALl-promoter [3l|, considering 
transcription activation by the TATA-binding protein. 
The operator position has been also shown to influence 
"bunchiness" of a promoter [33|. 

To the best of our knowledge, transcription bursts 
due to repressor as are analyzed in the present paper 
have not been observed yet experimentally. 

Transcription Interference: We have, at present, 
no direct experimental evidences for the possibility of 
de-repression by transcriptional interference(TI). Its bi- 
ological relevance, however, could be widespread in 
phage and E. coli; Convergent promoters is a common 
regulatory motif for all temperate phages that has a 
CII like protein, and about 100 examples of convergent 
promoters have been also found in E. coli [l3 | . 

To show how TI with de-repression could help us 
to understand a biological system, let us discuss the 
hyp-mutant of lambda; This system is intriguing be- 
cause of its high production of Cro in the lysogeny 
and its enhanced immunity against infection of other 
lambdas [bsI. Is^ . Its DNA configuration resembles that 
in Figdl^b) and the parameters in FigE] are matched 
to this system. Therefore, the maximal repression case 
there corresponds to the case of the promoter PR in 
lambda repressed by the factor of about 500 due to CI 
[ssl . [sgI . The strength of hyp-PRE is not known, 
but FigE] suggests that PR in lysogen could be de- 
repressed by the factor 10~30 due to TI from hyp-PRE, 
provided that open complex formation is fast and that 
CI binds to OR relatively slowly. This could explain, 
at least, a part of the large amount of Cro found in the 
lysogeny of the hyp- mutant. 

Another example is the Ofl30ij2 0ij3 mutant, which 
has been also found to show stable lysogens[38[, even 
though it is expected to be producing Cro 10^30 times 
more than a normal lambdaQ, as in the case of the 
hyp-mutant. Such similarity, i.e. the stable lysogens 
under the high production of Cro, between the On^iOpaO 
and the hyp-mutant leads us to speculate that the re- 
markable robustness of the lysogens [s^ of these phages 
should be rooted in the same unknown mechanism. 



Burst activity should be most directly monitored by a 
real time observation, but also can be examined quan- 
titatively from the number distribution of mRNA in a 
cell. This may be obtained if one can take snapshots of 
an assembly of cells from which the number of mRNA 
contained in each cell can be counted. The reaction 
rate constants for RNAP and TF should be able to be 
estimated from the mRNA distribution. 

For example, from the distribution one can calcu- 
late the Fano factor v, which is the ratio of the variance 
to the average 



{{n-{n)?) 
in) 



(39) 



with n being the number of mRNA in a cell. This 
Fano factor can be directly compared with our esti- 
mate of the number of transcriptions in a burst Tibst; 
In the bursting situation with nbst 3> 1, the Fano fac- 
tor should be given by ribst if the quiescent periods fol- 
low Poissonian process and are much longer than the 
bursting periods: 



«bst 



ko 



ko ~t~ ku 



for Tbst < TTF, (40) 



but is smaller than that if the burstings are not sepa- 
rated well enough: 



< "bst 



for Tbst < TTF- 



(41) 



In the case nbst <^ 1, we would have 1/ ~ 1 because 
each elongation initiation follows the Poissonian pro- 
cess. The full information of the distribution allows us 
more detailed comparison with our analysis. 

Another experiment we can propose is to construct 
DNA with a promoter exposed to a library of interfer- 
ing promoters with varying strength preferably in 
a parallel configuration to avoid RNAP collisions. Sup- 
pose the promoter is highly repressed by TF with un- 
known parameters. By examining how the promoter is 
de-repressed by the interfering promoters, the off-rate 
fcj^ of TF can be estimated as the lower limit of ^Ia 
that de-represses the promoter. 

Simplifications in the Present Treatment 

Before concluding, let us discuss some of the effects we 
have ignored in the present treatment. 



Roadblock: In our analysis of TI, we have assumed 
,hat RNAP always displaces TF without roadblock ef- 
'ect, but it is known that some TF's are roadblocks 
to RNAP. Roadblocks are most commonly reported in 
in- vitro experiments 4^, 41 , [4^ , whereas presence of 



elongation factors often allow RNAP to pass the road- 
block in in- vivo situations 4^, 4^. Reports on in- vivo 
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roadblocks is at present limited to the transcription fac- 
tor Lad and the restriction enzyme EcoRI 4J] . It has 
been reported that some roadblocks may be translo- 
cated, being pushed by two or more RNAP's 44l|. If 
two consecutive RNAP's are required to dislocate a 
TF, the activity of interfering promoter ^Ia in eq. ((37|) 
should be replaced by the effective activity, which is 
half of the original activity, 0^/2. This reduction fac- 
tor 1/2 should be further reduced in the case where 
a blocked RNAP may fall off before the second one 
arrives to give a push. 

Another possibility that RNAP does not remove 
TF is that RNAP just passes TF without displacing 
it; The repressor would simply not leave the vicinity of 
the operator, and thus maintains its function until it 
falls off by itself. This kind of situation has been actu- 
ally observed when an RNAP reads through a nucleo- 
some, displacing only parts of the histone complex [isl 
[jij . For some TF, one could imagine mixed situations, 
where TF is displaced but remains in physical proxim- 
ity during the RNAP passage. 

TP's such as CI in phage 186 ^ and CI in A on OR 
[stI do not act as roadblocks but are removed, which, 
we presume, are more common situations. 

Time difference for the promoter and the oper- 
ator: The interfering RNAP's clear/occlude the pro- 
moter pS first and then the operator in the convergent 
configuration, and the other way around in the parallel 
configuration. The time difference of the effects for the 
two sites depends on the distance between the two sites. 
If the binding times of TF or RNAP arc comparable 
or shorter than this time difference, we have to take 
this into account, which makes the situation favorable 
to the promoter(operator) in the convergent (parallel) 
configuration. 

RNAP collision: The RNAP from pS may be re- 
moved, even after it starts elongating, by colliding with 
the RNAP from pA. This effect is profound particularly 
when the distance between pS and pA is large. It has 
been found that the collision effect becomes substan- 
tial for convergent promoters with the pS-pA distance 
being of the order of v/ (2Ua), where 50 bp/sec) is 
the transcription elongation speed[3l- For the parallel 
configuration of promoters, the collision effect does not 
exist. 

Occlusion time: The occlusion time Tqcc, the time 
that the promoter pS is occluded by passing the RNAP 
fro m p A, was neglected. This has been also considered 
in [1J| and found that pS is influenced substantially by 
occlusion only when the activity of pA is stronger than 
0.1 



Mutual interference: In the convergent configura- 
tion of promoters, not only pA interferes with pS, but 
also pS interferes with pA. Such mutual interference 
effects are likely to be important in switching mech- 
anisms between equally strong convergent promoters, 
such as the convergent promoters PR and PRE of lambda 
phage in the early stages of infection. Full analytical 
treatment on the mutual interference is not easy in 
general case, but stochastic simulations [l3| and the 
four-world approximation analysis [l5| have been per- 
formed. In the present analysis, we consider the highly 
repressed promoter pS, thus the interference of pS on 
pA should be negligible. In the case of the parallel 
configuration, this effect does not exist. 
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Appendix: Promoter Activity in 
Michaelis-Menten Form 

Since the process of transcription initiation can be re- 
garded as an enzyme reaction, our results for the aver- 
aged promoter activity in the three step model can be 
put in the form of Michaelis-Menten kinetics. 

Let us start by the bare activity without TF. The 
binding rate ki, of RNAP should be proportional to the 
density of RNAP, 



kb = [RNAP] Kb 



(42) 



with a reaction constant k^,. Then, the bare activity 
([U can be written as 



[RNAP]/A-j^^AP omax 

^^°-l + [RNAP]/if*^^p''" 
with the maximum activity 

ko -\- fcg 

and the effective dissociation constant for RNAP 



K 



RNAP 



I, \ Qinax 
1 _|_ _ 1 



ko J l^h 



(43) 



(44) 



(45) 



TF has been introduced competitive inhibitor 
in our model. Its binding rate can be expressed as 



kr ^ mnr, 



(46) 
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with the TF density [TF] and the reaction constant 

lTF 



kF^ , then the dissociation constant for TF is given by 



TF 



(47) 



With these parameters, the expression ([25]) for the av- 
eraged activity with TF is written as 



TF 



[RNAP]/XS,NAP 



1 + [TF]/Atf + [RNAP]/A1nap 



n'l^'''', (48) 



which is in the standard form of Michaehs-Menten ki- 
netics with a competitive inhibitor. 
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Abstract 

Detailed derivations for the mathematical expressions in the text are given. 

I. Elongation Initiation Interval and Correlation 

Let C{t) be the time-dependent activity after the clearance of both the promoter and the operator. Then it 
can also be regarded as a correlation function of the transcription initiation, and it is related to the initiation 
interval distribution pij) as 



oo />oo 



C{t) = p{t) + dn dT2 S{t - n - T2)p(ti)p{t2) 

JO Jo 

/•oo /"OO /"OO 

+ / dri dT2 / drs (5(t - Ti - r2 - r3)p(ri)p(r2)p(r3) 
Jo JO JO 



(1) 



Since each term in the right hand side is a convolution of J'(t), the Laplace transformation 

C{s)= / C{t)e-''dt (2) 







can be obtained easily as a sum of geometrical series; 

p{s) 



with p{s) being the Laplace transform of p(t). 

II. Bare Promoters 

In this section, the explicit expressions for p{t) and C{t) for a bare promoter of each model are derived within 
the approximation that the self-occlusion effect is ignored. 

A. Single step model 

For the single step model, the elongation initiation is a simple Poissonian process with the rate fio, thus we 
have 

p{r) = n,e-''-\ p{s) = (4) 

and 

C{s) = ^, C{t) = no. (5) 
s 
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B. Two step model 

In the two step model, each elongation interval consists of an ofF-state period and an on-state period, whose 
length distributions, Poff (Toff) and Pon(Ton), are Poissonian given by 

Poff(Toff) = fcone~''°"^°", and Pon(Ton) = fcee~''"^°" , (6) 

respectively. Since the elongation interval is the sum of the off-period length and the on-period length, the 
elongation interval distribution p{t) for the two step model is given by 

/>oo pOC 

P{t) = / dToff / dTon(5(T - Toff - Ton)Poff(Toff)Pon(Ton)- (7) 

Jo Jo 

Again, the right hand side is a convolution of Poff ("Toff) and Pon('''on), thus in the Laplace transform, we have 



which gives 



Pir)={ fconfce j^^^^j^^ forfconT^fce _ 



fcg re for fcon = fee 

Then, the correlation function is given by 



= J(7Tf7TM' cw.«.(i-e-<--M.) (10) 

with Qq being the bare activity for the two step model: 

konke 1 _ 1 _ 1 , , 

"o^ fc +k ^ T +r ' = ^^^r- ^ ^ 

n^on ^ "-e ' on ^ ' e '^on "-e 

Note that the last expression simply represents that the average interval between elongation l/flo is the sum 
of the average waiting time to become the on-state Ton and the time for the elongation Tg. 



C. Three step model 

In the three step model, the initial transition between the off-state and the closed complex state is reversible, 
which means that the closed state goes either back to the off-state or forward to the open complex state with 
the branching ratios ku and kg, or with the probabilities 

^^VTJT' '''''^ '^ = VTir^^~^' ^^^^ 

respectively. Therefore, the promoter may get into the closed state many times before an RNAP starts elon- 
gation. Let n be the number of times that the promoter gets in the closed state before elongation, then the 
sequence of states and their probabilities are 

n state sequence probability 

(13) 



1 (off - closed) • open - elong q 

2 (off - closed) o (off - closed) • open - elong p q 



(off - closed o)" ^(off - closed) • open - elong p" ^q 



where o and • represent the branching probabilities p and q, respectively. 

Let Pn{T) be the elongation interval distribution for the interval during which the promoter goes through 
the closed state n times, then it is given by a convolution of the life time distribution of the off-state Pos{t), 
the closed state Pci{t), and the open state Pop(t); For example, Pi{t) is given by 

/•OO /"OO /"OC 

Pi{t) = dToff dTcl / drop (5 (r - Toff -Tcl -Top )Poff (Toff )Pcl(Tcl)Pop(7"op)- (14) 

Jo Jo Jo 
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In the same way, the Laplace transform of p„(t) for general n is given by 

Pn{s) = (poff{s)Pcl{s)j Pop{s) 



with 



Poff(s) 



Pop(s) 



(15) 



(16) 



s + kh s + ku + ko' s + kg 

The elongation interval distribution p{t) is the average over Pn{T) with the probability given by (jlSp . and 
can be calculated as follows; 



Pos{s)Pc\{s) 



is)pciis)p 



Popis)q = 



kffkgkg 



(s + fc+)(s + k_){s + ke) 



with 



* 2 



(fcfc + ku + ko) ± V (fch + /cu + ^o)^ - 4/cf,fco 



from which we obtain 
p(t) 



/lJ) 



(fee - fc+)(fce - A:_) 

From eqs.(I3]) and ([TT]), we have 



^{kb + ko + keY - ikbko 



.s{s + k^){s + k^)' 



kg 



{kb + ku + ko + fee) ± ^y {kb + ku + ko - keY 



4:khkr. 



which leads to 



cit) = no 



1 - 



k9e - k'^e 



fcVt 



where fio is the bare activity for the three step model: 
kbkokc 1 1 



kgk^ 



n 



n 



h* 



ko ko kb 



(17) 
(18) 
(19) 

(20) 
(21) 
(22) 



The average interval between elongations, l/r^o, are the sum of the three times: (1) r;,, the time for RNAP to 
bind and form the closed complex for the first time, (2) r*, the time for RNAP to form the open complex after 
the first binding, and (3) Tg, the time to start elongation. Note that the validity of this expression is not limited 
to the case within the two step approximation, where k*^ can be interpreted as the effective on-ratc. 

The time t* consists of two parts: (i) 1/fco, the time to go forward to the open state, and (ii) the re-binding 
time 1/fcb after unbinding multiplied by the average number of unbindings fc„/fco. Mathematical derivation of 
this expression is given in the appendix. We will encounter similar expressions in the following. 



III. Regulated Promoters 

Now, we derive the expressions for pir) and C{t) for each model of a bare promoter in the case where the 
promoter is repressed by a transcription factor (TF). In the case where the suppression is strong by a slow 
binding TF, the transcription activity occurs in bursts. The averaged activity is given by the long time limit 
t^oooi C[t). 

A. Single step model 

The state sequence between elongations can be classified according to the number of TF bindings, and the 
probability and the interval distribution for each case are obtained as 

n state sequence probability interval distribution 



off • elong q Poff(s) 

1 (oflF o TF) - off • elong pq Poff (s)pTF(s)Poff (s) 

2 (off o TF)2 - off . elong p'q {posis)PTFis)fpos{s), (23) 
n (off o TF)" - off • elong (poff(s)PTF(s))"poff (s) 
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where TF represents the state with TF at the operator site and 

L.T] 



Poff(s) 



fcTF + no 



k 

Ptf(s) = " 



l,TF 



no 



no 



Thus, we have 



Posis)q 



kr + no 



{s + kt')no 



n=0 



1 -Poff(s)pTF(s)p (.S + fc+)(s + fc-)' 



r TF , , TF 



and 



C{t) 



k^ + no± \/(fcTF + fcTF + noY - Ak^no 
no 



Z-TF I jLTF 



(fcTF ^ ^TFg-(fc-+fc-)t 



Thus, the steady state activity JItf for the single step model is given by 

i.TF 



f^TF = lim C{t) = 



no = 



?^bst 



Tbst + ttf 



with 



''bst = TtF ' "■bst = ^O'nDstj 



Ttf 



l,TF 



(24) 

(25) 
(26) 

(27) 

(28) 
(29) 



The last expression for fixF allows a simple interpretation in terms of bursting activity; Tbst and ttf are the 
bursting time and the quiescent time, respectively, and ribst is the number of transcriptions during the bursting 
time. 

B. Two step model 

For the two step model, the state sequences, their probabilities, and the interval distributions are 



n 




state sequence 




probability 


interval distribution 



1 

2 


(off c 
(off 


off • on - elong 
TF) - off • on - 
TF)2 - off • on - 


elong 
elong 


q 

pq 


Poff(s)Pon(s) 
PoS (s)PTF (s)Poff (s)Pon (s) 
(poff(s)PTF(s)) Poff(s)pon(s) 


n 


(off 


TF)" - off • on - 


- elong 


p^q 


(Poff (s)PTF (s)) "poff (s)pon (s) 



(30) 



with 

Poff(s) 



k'^^ "l~ ^on 
5 + ^ ~t~ ^on 



^TF ^ 



From these, we obtain 



7,TF 



kT^ ~t~ ko 



k^^ + fcon 



p{t) 



C{t) 



kakr, 



lTF 



(fc+ - ke){k- - ke) 



e ; ; — e + 



=^ 2 



TF , I TF 



fc-u + fcon i \/ (fc^ + fcon)^ ^ 4fcj|^F/j;^ 



TF 



Ai I ft/ -|- 



_ r,TF 



kl 



{k^^ + fcon + fc^^ + fce) i Y(fc^^ 4(A:onfcu^ + fcu^fce + fc^'^fce) 



(31) 

(32) 
(33) 
(34) 
35) 



The steady state activity JItf is now 



r^TF = lim C(t) 



"-bst 



fcon + fce + (fcb''^/fcu'^)fce 



Tbst + TT-F 



(36) 
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with 



"bst 



T'bst 



l,TF 



nhst- 



(37) 



The number of transcriptions during a bursting period is given by the winning ratio of RNAP to TF. 
It is interesting to see that eq. p6p can be also put in the form, 



r^TF — 



1 



lTF 



1 



1 



7'on(TF) = Ton + ^ ^TF, Ton = , Te — — , 



(38) 



which is similar to the expression for the bare promoter The expression for ron(TF) allows a similar 

interpretation with that for eqs. (f22|) : the time to reach the on-state from the off-state, Ton(TF) is the sum of 
the two times: (1) Ton, the time to reach the on-state without TF binding, and (2) the TF unbinding time, ttf, 
multiplied by the number of TF bindings, kj^ /kon, before an RNAP binds. 

In the case k"^^ ^ ke ^ k"^^ « 0, the correlation function C{t) shows a plateau. As the lowest order 
estimate, we put simply /cj^ ~ 0, then we obtain 



Cit) 



kpkr) 



lTF 



-k'rt 



k^ ~ fc, 



TF 



kr 



(39) 



therefore, C{t) behaves as 

Cit) « 

with the plateau time 
and the plateau value 



^on k^t 



jlTF 



for t < (A:TF + fcon)-i 
for (fcTF + fcon)"' <i <i 



pi 



kpkc 



jlTF 



■ exp 



' Z-TF I I. 



for t > t 



pi 



V = — (1 + 71bst) 



lTF 



• kp^ 



(40) 
(41) 



C. Three step model 

For the three step model, additional complication is that there arc two reversible transitions, i.e. the transition 
between the ofF-statc and the closed state, and the transition between the off-state and TF binding state, thus 
there exist two sequences of transitions within each elongation interval. This can be nicely represented by a 
binominal expansion; 

state sequence interval distribution with probability 

off cl • op - elong, 
( (ofr(X)TF-)-h(off0cl o) ) off cl • op - elong, 
( (off0TF-)-h(off0cl o) )%ff cl • op - elong, 

( (off0TF-)-h(ofr0cl o) )" off cl • op - elong. 



Poff(s)qiPci(s)g2Pop(s) 
(poff ( s )piPtf ( s ) + Poff ( s ) giPci ( S )P2 ) Poff ( s ) giPci ( s ) g2Pop ( s ) 

(poff(s)piPTF(s) +Poff(s)giPcl(s)p2) Poff(s)(7lPcl(s)g2Pop(s) 

(poff(s)piPTF(s) +Po«{s)qiPc\{s)p2Ypos{s)qiPc\{s)q2Pop{s) 



with the period length distributions 



k^ ~\~ kij ^ ^ ^ k^ _ ^ ^ k^ ~\- k, 

and the branching probabilities 



Poff(s) = , .TF , , ' PTf(s) = , " TF , Pc\{s) = —— Pop{s) = , 

s + k^'^ + kb s -f- fc,^ s + ku + ko s + ke 
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which are represented by the marks: ® for pi, for qi, o for p2, and • for q2. This gives 

■Poff{s)qiPci{s)q2Pop{s) 



p{s) = 



1 - (Poff(s)piPTF(s) + Pos{s)qiPci{s)p2) ' 



and the expression for p(t) is 



p{t) = kbkoke 



lTF 



(fci - ke){k2 - fce)(fc3 - ke) ^ (fc^+i - ki){ki^i ~ ki){ke - h) 

with —ki (i — 1,2,3) being the solution of the cubic equation 

+As'^ + Bs + C = 

with the coefficients 



A 
B 
C 



b + fcb + /c,jj + + ko 

TF/, , , , 7„ \ , 7„TF, 



kbkoky^ . 



Note that we define ki as a decay rate with a positive real part. 
From eqs.(|3]) and (|42|) . the correlation function C{t) is given by 



h.Ci.ChC 
'^1 '^2 '^3 



with — fcj- (i = 1, 2, 3) being the solution of the cubic equation 



with the coefficients 



D 

E 
F 



Ds^ + Es + F ^0 



— kf) ^" kb -{^ k^ -\- ky^ -\~ kQ -\- ke 

= {kj^ + + fco + ^e) + kb(k^^ + fco + ^'e) + {ku + ke,)ke 

= kbiko ~t" A^e) ~t~ (^b ~t~ ~t" ko^ke. 



The steady state activity I^tf is 

^TF = lun C(t) 

t — >-oo 

nbst 



fc^ A^o A^e k, 



TF 



k^^ kb{ko 4" A^e) ~t~ (^5^^ ~t" "1" ko)ke 



Thst + TTF 



with 



nhst 



ko 



kh ko 4^ 



Tbst 





for 


I.TF 

"■6 







?^bst 
ttf 


for 


i.TF 

w ~ 


-> OO 




1 








1 


Ttf + 

6 


"bst 


\,ko^ 


fee / ' 





(42) 



(43) 



(44) 



(45) 



(46) 



(47) 
(48) 

(49) 



and r^o being the bare activity of the three step model (j22p . The number of transcriptions ribst in a burst is now 
given by the winning ratio of RNAP kb/kj^ multiplied by the branching ratio in the closed state ko/{ko + ku). 
The expression for fixF can also be put in the form analogous to eq. ((22|) . 



TF 



1 

kb 



kl^ 



lTF 



1 

w 



- =T,(TF)+T:(TF)+Te (50) 



with 



y^TF ^ 

Tb(TF) = Tfc + • TTF, T:(TF) = To + • Tfc(TF), 

Kb r^o 



n 



kb ' 



(51) 
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This allows a similar interpretation with that for eq. (j38p ; The average interval of elongation 1 /fixF is the sum 
of three times: (1) ti,{TF), the time for RNAP to bind the promoter, (2) t*{TF), the time to form a open 
complex for the first time after RNAP binding, and (3) Tg, the time to elongate after forming the open complex. 
The time T^iTF) is the sum of (i) Tf,, the RNAP binding time, and (ii) the unbinding time, ttf, multiplied by 
the number of TF bindings, kj^ /kb, before an RNAP binds. Similarly, the time t*{TF) is the sum of (i) Tq, the 
time to form an open complex, and (ii) the binding time Tb{TF) multiplied by the average number of RNAP 
unbindings, ku/ko, before it forms an open complex. 



IV. Promoter Interference 

If there is another promoter, pA, competing with the promoter pS in the parallel or converging position, then 
their activities interfere with each other. Here, we consider only the interference effect on pS by pA. The activity 
of pA is given by the elongation interval distribution pA (t) . 

We consider two effects for the transcription interference: occlusion and sitting duck interference. The 
occlusion is the effect that an RNAP cannot bind to the promoter site of pS while the RNAP from pA is passing 
over the promoter site. The time that RNAP needs to pass through the promoter site is the occlusion time Tqcc- 
The sitting duck interference is that the RNAP sitting on the promoter site pS is removed by the RNAP from 
pA comes to pS. 

Under the interference of pA, the pS activity is limited within the elongation intervals from pA. Therefore, 
the average activity of pS under these effects fixi is the average over the activity within the interval r and the 
average by the probability that a time is in the interval of the length r; This is given by 

[ pa{t)t {- f °C{t) dt]dT 

SiTI = Too , (5^) 



PA{T)TdT 

using the transcription initiation correlation function C(t) without interference effects. 



A. Interference with unregulated promoters 

First, we consider the cases where the promoter pS is not regulated by TF. In this case, the correlation function 
C{t) is rather simple and the interference effect can be represented by a simple factor x, that is the averaged 
fraction of time that is not occluded: 



Pa{t)t ( —— ) dr 

(53) 



PA{T)TdT 

If we assume the simple Poissonian for pA with the activity VIa, 

Pa{t) = ^AC-^^^ (54) 

then X is given by 

XP = e-''^"-. (55) 

(i) In the case that pS can be described by the single step model, C(t) is given by the constant k^ = f^o as 
has been calculated ([5]), thus eq. (|5^ gives 

/ pa{t)t (- j f2o dt\ dr 

^Ti = ^ ^— = X^o, (56) 



/ PA{r)TdT 
Jo 



which means that the interference effects reduce the activity by the factor x- 

(ii) In the case of the two step model with the Poissonian pA, eq.([52| can be estimated as 
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using eq. pU]) . Here, fig is the bare activity ([TT|) for the two step modefl- 



B. Interference with regulated promoters 

Only difference from the cases above is that we use the correlation C{t) under the effect of TF. We give explicit 
expressions only for the Poissonian pA (|54[1 . 

For the single step promoter pS, using eas. (|28)) and (|52|) . wc obtain 



f^Ti — Xp ^0 



lTF 



l,TF 



i-TF I i-TF 



kr+kr 



XP ■ 



nhst 



Tbst 



(58) 



Ttf -r ^'■A 

The last expression simply shows that the quiescent 



with Xp given by eq. fSS]) and nbst and Tbst by eq, 
period is interrupted to be fi^^ when fi^^ < ttf- 

For the two step promoter, we give the expression only for Tqcc = 0, namely, without the occlusion effect 



ke kon k^ 



TF 



keko 



k^-k^ 



flA 



lTF 



^A 



k^ fiA + k^^ k^l + 



TF 



^Ai^A + k% + kC- k^k^/kY) 

i^A + k^){nA + k^) 



with given by eg. ([55]) and JItf being the steady activity under TF by cq. ipSl) . 

In the same approximation as cq. (|39[) . kj^ S> fee ^ fcjf^ ~ 0, using the approximate form of k^, we obtain 



n 



nhst 



TI 



fcjj^^ + ko 



^A + kj^kjik'^^ + fcon) 1 + (fcTF + k,,,)/n, 



rihst^A 

»bst 

tpl 

Ibst 
tpl 



for Ha <. l/tpi 



TF 



^hst kT^ + ko 



^A 



h h 



for 1/tpi < < (fcfc" + fcon) 

for (fcTF + kon) < 



with ribst (|37p and tpi (j40p . The maximum value of ilxi is achieved at VLa ~ a/ n.^ 
([39|) and JIti as a function of correspond to each other by the correspondence t ^ 1 /VI a 
For the three step model, we give only a formal solution for Poissonian pA: 

3 



kj^ke. The behavior of C{t) in 



h,hh i-TF 
"TI = , r',n k^koke 



i^^-k^) 



VLa 



"■1 "-2 '^3 



t A:f (A:f+, - fcf )(fcf_, - kf) + fcf ' 



(59) 



with the same kf {i = 1, 2, 3) with those in eg. (1^5]) . 



Appendix: Mathematical explanation for r* in (1221) 



In this appendix, we will give a mathematical explanation for the expression of r* in cq. (j22p . 

This is the time for RNAP and promoter to form the open complex after the first binding of RNAP. Since 
the initial binding/unbinding process is reversible, after the first binding, the system may either go forward to 
form the open complex, or may go backward to unbind with the probabilities. 



ko ^" ku 



and p 



ko ~t~ ku 



(A.l) 



^ Ea. l|57| l disagrees with the corresponding expression 



X(fcon + fce) 



of eq.(4) and Figure 1(d) in Sneppen et al. [J. Mol. Biol. 346 (2005) 399— 409] (the notations have been changed from the original 
ones). It is not difficult, however, to see that this expression cannot be correct because this does not reduce to that for the single 
step case I I56I I in the fcon — > oo limit. In its derivation, only the on-rate was reduced by the factor x, and it was not taken into 
account that the total probability of the states available to RNAP is limited by the factor x- 
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respectively. The average waiting time for either of the cases to happen is 

1 



(A.2) 



If the system goes backward to unbind the RNAP, then after the thne 

another RNAP binds to form the closed complex again, and the situation becomes the same as before. 

If the system goes backward to unbind n times before it proceeds to form the open complex, the times that 
the system spends and the probabilities that should occur arc given 

n time probability 

Ivq -\- kij /lq -\- (A. 4) 



1 ^1 

ko -\- ku kjy 



p q 



thus the average time is given by 



oo 

E 

n=0 



1 



ko H~ ku 



1 

kb 



ko 



g — — + — — 

ko ko kb 



(A.5) 



which is T* in eq. 



